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ABSTRACT 


Projectile  Energy  Dependence  of  Aluminum  and 
Silicon  Ka  X-Ray  Satellites.  (August  1978) 
Blake  Isamu  Sonobe,  B.  S.,  USAF  Academy 
Chairman  of  Advisory  Committee:  Dr.  Rand  L.  Watson 

7 K ftyM 


The (Kq  satellite  spectra  of  aluminum  and  silicon  produced  by 
collisions  with  5.4  to  40.6  MeV  helium  ions,  12.0  to  74.1  MeV  carbon 
ions,  and  23.4  to  117.0  MeV  neon  ions  have  been  measured.  The  apparent 
average  L-vacancy  fractions  <(p  ) determined  from  the  relative 


intensities  of  the  satellite  peaks  were  observed  to  be  dependent  on  the 
projectile  energy  and  atomic  number  and  on  the  target  atomic  number. 
Comparisons  of  the  binary-encounter  approximation  (BEA)  and  semi- 
classical  approximation  (SCA)  predictions  with  the  experimental  p 

JL 

values  as  a function  of  the  projectile-to-L-shell  electron  velocity 
ratio  for  helium  on  aluminum  are  presented.  It  was  found  that  the 
BEA  greatly  over-predicted  the  degree  of  L-shell  ionization  and  failed 
to  accurately  predict  the  shape  of  the  velocity  dependence.  The  SCA 
predicted  the  shape  of  the  dependence  at  the  higher  velocities  but  was 
too  low  in  magnitude. 
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CHAPTER  I 
INTRODUCTION 

Since  the  discovery  of  x-rays  in  1895  by  Roentgen,1  who  observed 

this  mysterious  radiation  from  radioactive  sources,  many  investigations 

have  been  made  to  produce  and  characterize  x-rays.  In  the  early  1900's 

x-rays  were  produced  by  bombardment  of  various  targets  with  alpha 

2-i2 

particles  from  radioactive  sources  and  by  electrons.  In  the  late  ^ 

13 

1920 's,  Gerthsen  performed  experiments  using  protons  and  alpha 

14 

particles  to  study  the  production  of  inner-shell  vacancies.  Coates, 
in  1934,  found  x-rays  resulting  from  heavy  ion  bombardment  to  be 
characteristic  of  the  atoms  from  which  they  were  emitted.  One  of  the 
early  studies  of  x-rays  produced  by  charged  particles  from  an  accelera- 
tor came  in  1953  when  Lewis1'*  observed  tremendous  amounts  of  character- 
istic x-rays  from  proton  bombardment  (energies  ranging  from  1 to  3 MeV) 
on  thick  tantalum  targets.  Huus  and  Zupancic1^  also  observed  large 
amounts  of  characteristic  x-rays  while  studying  gamma  rays  from  proton 

bombardment.  , 

Energetic  charged  particles  from  accelerators  have  since  been 
frequently  used  to  produce  x-rays.  The  passage  of  a heavy  charged 
particle  through  matter  often  results  in  the  ejection  of  one  or  more 
of  the  inner-shell  electrons  from  the  atomic  orbitals  of  the  target 
atom.  This  ionization  is  characterized  by  the  emission  of  x-rays  or 


The  citations  in  this  thesis  follow  the  style  of  the  Physical 
Review  A. 
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Auger  electrons  during  the  subsequent  decay  of  the  excited  atom.  Ka 

x-rays  are  emitted  as  a result  of  a 2p->-ls  electron  transition  in  an  atom 

with  one  K-shell  vacancy.  Ka  x-ray  satellites  result  when  this 

transition  occurs  in  an  atom  having  one  or  more  L-shell  vacancies  in 

addition  to  the  K-shell  vacancy  at  the  time  of  x-ray  emission.  The 

absence  of  one  or  more  L-shell  electrons  reduces  the  screening  of  the 

nuclear  potential  experienced  by  the  remaining  electrons  and  increases 

their  binding  energies.  The  Ka  x-ray  transition  energy  is  increased 

above  its  characteristic  or  "normal"  value  and  produces  Ka  x-ray 

satellite  peaks  in  the  x-ray  spectrum. 

9 

In  1916,  Siegbahn  and  Stenstrom  first  noticed  x-ray  satellites  as 
low  intensity  lines  at  slightly  higher  energies  than  the  normal  Ka  „ 

y 

line.  Their  measurements  extended  from  sodium  (atomic  number  11)  to 

17  18 

zinc  (atomic  number  30).  In  1921,  Wentzel  * postulated  that  these 

lines  originated  from  atoms  having  one  or  more  L-shell  vacancies  in 

19 

addition  to  the  K-shell  vacancy.  Richtmyer  proposed  that  the 
satellite  lines  were  the  result  of  two  simultaneous  electron  jumps 
resulting  In  a single  x-ray.  The  Wentzel  theory  has  since  been  proven 
to  be  the  correct  one. 

20*~  2 A 

More  recently,  many  authors  have  studied  multiple  inner-shell 

ionization  by  heavy  ion  projectiles  by  observing  the  shifts  of  the  Ka 
and  KB  x-ray  energies.  They  noted  that  as  the  degree  of  L-shell 

ionization  increased,  the  energy  shifts  became  correspondingly  larger. 

/ 25-28 

High  resolution  Bragg  spectrometers  have  enabled  investigators  to 

resolve  the  Ka  x-ray  satellite  lines.  Various  theoretical  calculations 

29  30  31 

(Herman-Skillman  HFS , C.  Froese  Fischer  HF,  Veigele  program,  etc.) 


T 
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27 

have  confirmed  the  identity  of  these  linr  . McCrary  and  Richard  also 

observed  multiplet  structure  in  the  Ka  x-ray  spectra  of  silicon  produced 

32 

by  bonfcardment  with  30  MeV  oxygen  ions.  As  early  as  1929,  Ray, 

33  34 

Langer,  and  Wolfe,  using  spectroscopic  interpretations,  predicted 

that  satellites  were  composed  of  a number  of  multiplet  lines. 

35  36  37 

Bhalla  * and  Kauffman  et  al.  have  attempted  to  assign  relative 
intensities  to  the  various  spectroscopic  terms.  An  illustration  of  the 
complexity  of  the  Ka  x-ray  satellite  spectrum  of  neon  is  shown  in 
Fig.  1.  The  spectra  in  Fig.  1 were  constructed  by  assuming  a statistical 
population  of  the  spectroscopic  terms  and  a binomial  distribution  for 
the  relative  population  of  the  electronic  configurations  (p  = 0.375), 

L 

and  by  using  the  x-ray  transition  energies  and  multiplet  partial 

35 

fluorescence  yields  of  Bhalla.  The  relative  intensities  of  the 
multiplet  components  for  the  various  2s-  and  2p-orbital  (L-shell)  vacancy 
configurations  are  given  in  the  lower  three  spectra.  The  top  spectrum 
is  a composite  of  these  three  spectra.  The  notation  (KLn)  used  in  the 
top  spectrum  denotes  configurations  with  one  K-shell  and  n L-shell 
vacancies.  Despite  the  overlapping  multiplet  components,  most  lines 
are  grouped  together  such  that  individual  peaks  result.  These  multiplet 
components  and  the  reduced  screening  of  the  nuclear  potential  resulting 
from  M-shell  vacancies  in  addition  to  the  K-  and  L-shell  vacancies 
existing  at  the  time  of  x-ray  emission  contribute  to  the  broadening  of 
the  peaks. 

Typical  silicon  Ka  x-ray  spectra  produced  by  2 MeV/amu  He,  C,  0,  Ne, 
S,  and  Ar  ions  are  shown  in  Fig.  2.  The  first  peak  on  the  left 
corresponds  to  the  normal  Ka  9 x-ray  line.  Each  succeeding  peak  to  the 


RELATIVE  COUNTS /CHANNEL 


ENERGY  (keV) 

FIG.  2.  Silicon  Ka  x-ray  satellite  spectra  produced  by 
2 MeV/amu  ions. 
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right  is  a satellite  line  associated  with  Ka  x-ray  emission  in  the 
presence  of  one  through  six  L-shell  vacancies.  The  KL^  peak  is  cut  off 
in  the  S and  Ar  spectra  by  the  silicon  K absorption  edge  which  occurs 
at  1838  eV . The  intensity  distribution  rapidly  shifts  to  the  higher 
order  satellites  as  the  atomic  number  of  the  projectile  increases.  The 
relative  intensities  of  the  satellites  in  the  He,  C,  0,  and  Ne  spectra 

OO 

can  be  well  described  by  binomial  distributions.  This  appears  to  be  a 

characteristic  of  satellite  spectra  produced  by  ions  having  atomic 

39  40 

numbers  less  than  that  of  the  target.  ’ Departures  from  binomial 

distributions  caused  by  the  enhancement  of  the  Ka  peak  are  observed 

1 » *■ 

in  the  S and  Ar  spectra.  This  feature  is  only  observed  when  the 

projectile  atomic  number  is  greater  than  the  target  atomic  number. 

The  mechanisms  contributing  to  this  enhancement  are  not  fully  understood. 

Detailed  studies  have  shown  that  the  Ka  x-ray  satellite  structure 

42  43 

is  dependent  upon  the  charge  state  of  the  projectile,  ’ the  chemical 

39  44-47 

environment  of  the  target  atom,  ’ the  atomic  number  of  the 

projectile, the  atomic  number  of  the  target  atom,^ 
and  the  energy  of  the  projectile. »54  56  the  atomic  number  of 
the  projectile  is  increased,  the  degree  of  L-shell  ionization 
increases,  and  as  the  atomic  number  of  the  target  atom  is 

increased,  the  degree  of  L-shell  ionization  decreases . ^ 

Kauffman  et  al.  ’ have  observed  in  studies  of  carbon,  nitrogen, 
oxygen,  and  fluorine  ions  on  neon  that  as  the  charge  state  of  the 


projectile  increases,  the  degree  of  multiple  L-shell  ionization 
increases.  In  the  chemical  effects  investigations,  energy  shifts 


44-46 


of  the  Ka  and  Kf5  structures  and  shifts  in  the  relative  Intensities  of 
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39  47 

the  satellites  ’ were  noted. 

Despite  the  numerous  studies  of  Ka  satellite  dependence  on  the 

projectile  energy  carried  out  in  recent  years,  a comprehensive 

investigation  encompassing  a wide  energy  range  and  compared  with  the 

predictions  of  an  improved  binary-encounter  approximation  (BEA)"^ 

58 

formulation  and  with  the  semi-classical  approximation  (SCA)  has  not 
been  accomplished.  The  applicability  of  the  SCA  to  the  prediction  of 
single  K-  and  multiple  L-shell  ionization  probabilities  has  not  been 
tested.  Comparison  of  experimental  data  with  the  predictions  of  the 
BEA~^  was  done  by  Olsen  and  Moore. They  used  24.0  to  48.0  MeV  oxygen 
ions  to  study  the  calcium  Ka  x-ray  satellite  and  hypersatellite 
structures.  The  BEA  predicted  the  shape  of  the  ionization  cross  section 
dependence  on  the  projectile  energy  but  was  approximately  50  percent 

54 

higher  than  the  experimental  values.  In  other  studies,  Knudson  e£  al. 

measured  the  aluminum  Ka  x-ray  yields  from  collisions  with  neon  ions 

over  an  energy  range  from  1.5  to  15  MeV,  however,  no  comparisons  were 

made  with  theoretical  models.  The  degree  of  L-shell  ionization  in  these 

49 

measurements  increased  with  increasing  projectile  energy.  Li  «Jt  al . 

noted  that  the  fraction  of  KB  x-rays  emitted  in  the  presence  of  an 

L-shell  vacancy  decreased  with  increasing  projectile  velocity  above  a 

projectile-to-L-shell  electron  velocity  ratio  of  1 for  deuterons  and 

alpha  particles.  This  study  was  the  first  to  test  the  applicability  ?f 
59  57 

the  Gryzinski  and  BEA  models  to  the  description  of  single  K-  plus 
single  L-shell  ionization.  Richard  elt  al.  ^ measured  the  aluminum  Ka 


x-ray  yields  using  0.4  to  3.0  MeV  helium  ions.  They  compared  the 
measured  ionization  cross  sections  with  the  predictions  of  a modified 
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BEA  model  and  noted  that  the  gross  features  of  the  data  were  fairly 

well  explained  but  several  areas  of  disagreement  existed.  This  model 

assumed  that  the  single  K-  and  multiple  L-shell  ionization  cross  section 

52 

is  given  by  a binomial  distribution.  Most  recently.  Hill  et  al.  did 
a comprehensive  study  of  the  projectile  atomic  number  and  energy 
dependence  of  titanium  Ka  x-rays  and  satellites.  Beams  of  2.5  to  5.0 
MeV  protons,  4 to  18  MeV  helium  ions,  10.5  to  26.25  MeV  lithium  ions, 

21  to  45  MeV  carbon  ions,  and  15  to  60  MeV  oxygen  ions  were  used  to 
excite  thick  titanium  targets.  With  the  exception  of  the  oxygen  ion 
data,  the  relative  intensities  of  the  satellites  produced  by  the 
projectiles  used  in  this  study  showed  rough  qualitative  agreement  with 
the  predictions  of  a theoretical  model  by  Hill  and  Madison  (discussed 
in  Ref.  52).  There  was  poor  agreement  between  the  oxygen  ion  data  and 
the  theoretical  prediction.  This  model  also  assumed  a binomial 
distribution  for  the  multiple  ionization  cross  sections  and  used 
empirical  correction  factors  to  compensate  for  underpredictions  at  small 
impact  parameters. 

The  objective  of  this  study  then  is  to  test  the  applicability  of 
the  SCA  and  BEA  to  predict  single  K-  and  multiple  L-shell  ionization 
over  a wide  energy  range.  An  extension  of  the  investigation  by  Richard 
et  al. ^ (0.4  to  3.0  MeV  He  on  aluminum)  to  much  higher  energies  was 
made.  The  aluminum  and  silicon  Ka  x-ray  spectra  resulting  from  impacts 
with  5.4  to  40.6  MeV  helium  ions  were  measured  and  the  aluminum  results 
compared  with  the  predictions  of  the  BEA  and  SCA.  The  Ka  x-ray  satellite 
dependence  on  the  energy  of  heavy  ion  projectiles  was  also  studied  using 
21.0  to  74.1  MeV  carbon  and  23.4  to  117.0  MeV  neon  ions. 
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In  Chapter  II,  the  projectiles,  targets,  and  experimental  procedures 
used  to  produce  multiple  inner-shell  vacancies  and  to  detect  the  x-rays 
are  described.  The  method  of  data  analysis  is  presented  in  Chapter  III 
and  the  results  of  this  study  given  in  Chapter  IV.  Finally,  three 
prevailing  theories  of  inner-shell  ionization,  the  plane-wave  Bom 
approximation,  the  semi-classical  approximation,  and  the  binary- 
encounter  approximation,  will  be  discussed  briefly  in  Chapter  V and  the 
results  of  the  SCA  and  BEA  compared  with  experimental  data.  The 
Appendices  contain  a computer  listing  for  the  Gaussian  integration 
routine  used  to  calculate  the  theoretical  ionization  cross  sections  for 
the  various  vacancy  configurations  and  a discussion  of  the  dependence  of 
the  sensitivity  of  x-ray  detection  on  the  atomic  number  and  energy  of 
the  projectile. 
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CHAPTER  II 

EXPERIMENTAL  PROCEDURES 
A.  Beams  and  Targets 

The  Texas  A&M  University  variable  energy  cyclotron  was  used  to 
accelerate  helium  ions  to  energies  of  5.4,  6.6,  10.4,  23.2,  and  40.6 
MeV;  carbon  ions  to  energies  of  12.0,  22.0,  26.9,  42.9,  and  74.1  MeV; 
and  neon  ions  to  energies  of  23.4,  41.6,  and  117.0  MeV.  The  beams  were 
extracted  from  the  88-inch  cyclotron  and  directed  down  an  evacuated  beam 
line  to  the  experimental  area  by  two  switching  magnets.  They  were 
aligned  and  focused  to  a spot  approximately  0.6  cm  in  diameter  using  two 
sets  of  quadrupole  magnets.  Alignment  and  focusing  was  observed 
visually  by  allowing  the  beam  to  strike  a fluorescent  material  (a  mix- 
ture of  zinc  sulfide,  ZnS,  and  cadmium  sulfide,  CdS)  placed  in  the 
target  position  and  adjusting  the  current  on  the  quadrupole  magnets 
while  monitoring  the  beam  spot  with  a closed  circuit  television  camera. 
The  camera  was  located  at  the  side  of  the  target  chamber  to  view  the 
fluorescent  target  through  a plexiglas  panel  on  the  target  chamber.  The 
camera  could  be  controlled  remotely  from  the  cyclotron  control  room. 
After  the  beam  had  been  aligned  to  the  center  of  the  target  area  and 
focused,  a 0.6  cm  diameter  collimator  was  remotely  positioned  at  the 
front  of  the  target  chamber  and  the  fluorescent  target  replaced  by  a 
target  of  experimental  Interest.  The  target  assembly  consisted  of  a 
target  wheel  capable  of  holding  three  targets  plus  the  fluorescent 
target.  The  target  wheel  was  remotely  controlled  and  a target  was 
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selected  by  simply  rotating  the  target  wheel  until  the  desired  target 
was  in  the  proper  target  position.  The  beam  current  measured  in  the 

i 

target  was  maintained  between  50  and  100  nanoamps . 

The  aluminum  and  silicon  Ka  spectral  measurements  were  performed 

using  thin  aluminum  and  silicon  monoxide  and  thick  silicon,  silicon 

monoxide,  and  silicon  dioxide  targets.  The  aluminum  targets  consisted 

2 

of  self-supporting  aluminum  foils  approximately  2.34  mg/cm  thick  glued 

to  2.54  cm  diameter  aluminum  rings.  The  thin  silicon  monoxide  targets 

were  prepared  by  vacuum  evaporation  of  silicon  monoxide  onto  thin 

aluminum  backings,  and  the  backings  were  then  glued  to  2.54  cm  diameter 

aluminum  rings.  The  thicknesses  of  these  silicon  monoxide  targets 

2 2 

ranged  from  0.92  mg/cm  to  1.71  mg/cm  . The  thicknesses  given  above  for 
the  aluminum  and  silicon  monoxide  targets  are  the  effective  thicknesses 
and  have  been  corrected  to  account  for  the  45°  inclination  of  the  target 
with  respect  to  the  beam.  The  thick  silicon  target  was  a silicon 
crystal  taken  from  a retired  Si (Li)  detector  and  mounted  on  a 2.54  cm 
diameter  aluminum  ring.  The  thick  silicon  monoxide  and  silicon  dioxide 
targets  were  prepared  by  compacting  finely  ground  powders  in  a pellet 
press,  the  resulting  pellets  (1.3  cm  diameter,  0.3  cm  thick)  were  glued 
to  aluminum  backings,  and  the  backings  were  then  glued  to  2.54  cm 
diameter  aluminum  rings. 

B.  X-ray  Measurements 

A high  resolution  plane  crystal  spectrometer  was  used  in  these 
experiments  to  obtain  the  Ka  satellite  spectra  of  aluminum  and  silicon. 

The  spectrometer  had  been  adapted  for  in-beam  use  with  a high-energy 
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particle  accelerator  and  was  constructed  from  a design  obtained  from 

the  Naval  Research  Laboratory.^1  It  Is  essentially  Identical  to  the 

45  25 

spectrometer  used  by  Burkhalter  e£  al.  and  Knudson  eit  al.  in 
previous  measurements  of  ion-induced  x-ray  spectra. 

As  shown  in  Figs.  3 and  4,  the  crystal  spectrometer  consists  of 
a flat  crystal  mounted  at  the  center  of  a 15.24  cm  turntable  on  which 
the  detector  rides.  A relationship  is  maintained  between  the  crystal 
mount  and  the  turntable  such  that  as  the  mount  rotates  through  the  Bragg 
angle  0 , the  turntable  rotates  through  the  angle  20 . A gear  train 
driven  by  a percision  stepping  motor  maintains  the  relationship  between 
the  crystal  mount  and  the  turntable.  The  detector  is  thus  in  the 
correct  position  to  observe  x-rays  which  satisfy  the  Bragg  relationship 

nA  - 2d  sin  0,  (2.1) 

where  2d  is  double  the  atomic  spacing  of  the  crystal,  A is  the  wave- 
length of  the  x-ray,  0 is  the  angle  of  reflection  from  the  crystal,  and 
n is  the  reflection  order.  Modifications  to  the  spectrometer  to  allow 
adjustment  of  the  Bragg  crystal  and  the  entrance  soller  slit  alignment 
remotely  can  be  seen  in  Fig.  4. 

The  spectrometer  was  positioned  perpendicular  to  the  incident  ion 
beam,  and  the  target  was  mounted  at  45°  with  respect  to  both  the  ion 
beam  and  the  spectrometer.  The  incident  x-rays  were  collimated  prior 
to  impinging  on  the  Bragg  crystal  by  a set  of  commercially  available 

soller  slits  having  a length  of  3.18  cm  and  a plate  spacing  of  0.127  mm 

which  limited  the  angular  divergence  of  the  x-rays  to  approximately 
0-3°.  A similar  set  of  soller  slits  having  a plate  spacing  of  0.254  mm 

was  used  in  front  of  the  detector  to  reduce  the  background  due  to 


FIG.  3.  Crystal  spectrometer  experimental  configuration. 
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FIG.  4.  Crystal  spectrometer  with  crystal  and  soller  slit 
adjustments. 
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scattered  x-rays  and  electrons.  Since  the  best  resolution  is  obtained 
with  a crystal  having  a 2d  spacing  only  slightly  larger  than  the  wave- 
length of  interest,  a flat  PET  (pentaerythritol)  crystal  was  selected 

o 

for  these  measurements.  The  2d  spacing  for  a PET  crystal  is  8.75  A 

o 

and  7.13  A respectively.  A gas  flow  proportional  counter  (90%  argon 

2 

and  10%  methane)  having  an  aluminum  coated  (100  ug/cm  ) Mylar 
2 

(530  pg/cm  ) window  was  used  to  detect  the  reflected  x-rays.  Lead 
shielding  was  placed  between  the  collimators  and  the  spectrometer  to 
reduce  the  background  from  particles  scattered  by  the  collimators. 

Automatic  control  of  the  spectrometer  was  provided  by  a specially 
designed  system.  A schematic  diagram  of  the  experimental  layout  is 
shown  in  Fig.  5.  Signals  from  the  -.-itector  were  amplified  by  a 
Canberra  Model  1417  spectroscopy  amplifier  and  sent  to  a Canberra  Model 
1437  timing  single  channel  analyzer  where  an  energy  window  was  set 
around  the  x-rays  of  interest.  The  output  signals  from  the  single 
channel  analyzer  were  then  sent  to  the  spectrometer  control  where  an 
internal  gate  permitted  passage  of  the  signals  to  a Northern  NS-6300 
multichannel  pulse  height  analyzer  operated  in  the  external-multiscale 
mode  and  to  a PDP  15  computer.  Passage  of  the  signals  was  restricted 
each  time  the  spectrometer  rotated  to  a new  counting  position  and  while 
the  multiscalers  on  the  pulse  height  analyzer  and  the  PDP  15  were  being 
advanced  and  reset. 

A block  diagram  of  the  spectrometer  control  system  is  shown  in 
Fig.  6.  The  counting  time  at  each  spectrometer  angle  was  determined  by 
a r reset  number  of  monitor  counts  selected  on  thumbwheel  switches 


located  on  the  front  of  the  spectrometer  control.  The  number  of  monitor 


) 


FIG.  5.  Crystal  spectrometer  experimental  layout. 
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counts  was  selected  such  that  approximately  3000  counts  above  the  back- 
ground were  obtained  in  the  peak  channel  of  the  satellite  of  greatest 
intensity.  The  monitor  counts  were  generated  by  a digital  current 
integrator  which  integrated  the  current  produced  in  the  target  by  the 
beam.  When  the  selected  number  of  monitor  counts  had  been  collected, 
the  spectrometer  control  system  automatically  stepped  the  spectrometer 
to  the  next  counting  position  and  advanced  and  reset  the  multiscalers 
on  the  pulse  height  analyzer  and  the  PDP  15  computer.  The  number  of 
steps  between  one  counting  position  and  the  next  can  also  be  selected 
by  another  set  of  thumbwheel  switches  on  the  front  panel  of  the 

_3 

spectrometer  control  system  (9.09  x 10  degree/step).  Both  the  number 
of  monitor  counts  being  accumulated  and  the  spectrometer  counting 
position  were  displayed  on  the  front  panel  of  the  spectrometer  control 
system  by  arrays  of  light  emitting  diodes. 

The  PDP  15  computer  provided  the  capability  of  on-line  data 

i 

analysis  of  a previously  stored  spectrum  while  a new  spectrum  was 
accumulating.  The  data  analysis  program  was  adapted  for  the  PDP  15 
L during  the  progress  of  these  experiments  and  the  PDP  15  was  used  in 

| the  later  measurements. 

I 
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CHAPTER  III 
DATA  ANALYSIS 


The  x-ray  signals  coming  from  the  spectrometer  control  system  to 
the  pulse  height  analyzer  were  recorded  in  one  of  256  channels  of  the 
analyzer  memory  or  into  one  of  the  storage  locations  of  the  PDP  15 
computer.  Each  channel  or  storage  location  corresponded  to  one  of  the 
spectrometer  counting  positions.  Prior  to  the  use  of  the  PDP  15,  once 
a spectrum  was  completed,  it  was  copied  onto  a magnetic  tape  by  means  of 
a Kennedy  Incremental  1600  tape  drive  unit  which  was  used  in  conjunction 
with  the  pulse  height  analyzer.  The  tape  was  processed  on  an  IBM  7094 
computer  using  a program  named  DIABOLICAL  PLOT  which  had  the  capability 
of  reading  data  tapes  from  several  analyzer  systems,  producing  punched 
cards  containing  spectral  data,  and  graphically  representing  each 
spectrum  by  printer  plots  and/or  calcomp  plots  (linear  or  semi-log). 
Using  the  punched  data  cards,  each  spectrum  was  further  analyzed  by 
carrying  out  a non-linear  least-squares  fit  to  the  spectrum  employing 
a Gaussian-plus-exponential-tail  peak  fitting  program  called  AIRLIFT  to 
obtain  the  relative  satellite  intensities  and  peak  positions.  Further 
details  and  discussions  of  this  peak  fitting  procedure  can  be  found 
elsewhere.  An  example  of  the  fits  typically  obtained  is  shown  in 

Fig.  7.  The  solid  circles  represent  data  points  and  the  open  circles 
represent  calculated  points. 

The  first  peak,  the  normal  Ka1  _ peak,  was  assigned  the  Ka..  - 
energy  given  by  Bearden. ^ This  energy,  ER,  is  related  to  the  Bragg 
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FIG.  7.  Example  of  a typically  fitted  spectrum. 


angle,  0 , by  the  expression 
i\ 
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(3.1) 


where  n = 1 for  first-order  reflections  and  where  = — has  been 

* h 

substituted  into  Eq.  (2.1).  In  this  substitution,  is  the  wavelength 

of  the  Ka  x-ray,  h is  the  Planck  constant,  and  c is  the  speed  of  light. 

Any  succeeding  position  on  the  spectrum  corresponds  to  0=0-  A0 

IK.  1 

where  A0^  can  be  determined  by  knowing  the  number  of  channels  from  the 
Ka1  „ position,  the  number  of  steps  per  channel  selected  on  the 
spectrometer  control,  and  the  number  of  degrees  per  step  of  the 
spectrometer.  The  corresponding  energy  is  then 


nhc 

2d  sin  0, 


(3.2) 


Each  peak  intensity  was  corrected  for  absorption  in  the  target, 
absorption  in  the  proportional  counter  window,  and  for  detection 
efficiency  according  to  the  equation 


T T T 
*12  3 


(3.3) 


where  I is  the  corrected  peak  intensity,  is  the  number  of  x-rays 
detected  in  the  photopeak,  is  the  correction  for  x-ray  absorption  in 
the  target,  is  the  correction  for  x-ray  absorption  in  the  detector 

window,  and  is  the  detector  efficiency  for  the  photopeak.  The  terms. 


T^  and  , were  determined  from  the  expression 
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where  u is  the  mass  absorption  coefficient  of  the  target  or  the  detector 

66 

window  taken  from  a compilation  by  Storm  and  Isreal  and  x is  the 
effective  target  or  detector  window  thickness.  The  detection  efficiency 
was  determined  by 

T3  « 1 - e_ux  , (3.5) 

where  u is  the  mass  absorption  coefficient  of  the  detector  gas  and  x is 
the  effective  thickness  of  the  detector  gas.  Since  relative  intensities 
and  not  absolute  intensities  were  measured,  corrections  for  crystal 
reflectivity  efficiency  were  not  made.  Also,  the  crystal  reflectivity 
is  not  expected  to  change  significantly  over  the  range  of  angles  being 
scanned  for  each  spectrum  (65.4°  to  73.2°  for  aluminum  and  46.9°  to 
50.3°  for  silicon).  The  corrections  for  absorption  by  the  target  and 
detector  window  and  for  detection  efficiency  had  a very  small  effect 
since  only  relative  intensities  were  determined  and  the  energy  span  of 
each  spectrum  was  less  than  120  eV. 

The  effective  thickness  for  each  projectile-target  confoination  was 
determined  by  calculating  the  average  depth  of  x-ray  emission  from  the 
front  of  the  target.  The  rate  per  unit  solid  angle  of  x-ray  emission 
from  the  front  of  the  target  is 
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(3.6) 


where  is  the  incident  particle  rate,  co  is  the  fluorescence  yield,  f 


is  the  weight  fraction  of  the  element  of  interest,  N is  Avogadro's 

A 

number,  W is  the  atomic  weight  of  the  element  of  interest,  t is  the 

al  • 

2 

target  thickness  in  g/cm  , a(x)  is  the  cross  section  for  K-shell 
ionization  at  depth  x,  and  T(x)  is  the  transmission  correction  given 

— V1X 

by  T(x)  = e . The  values  of  o(x)  were  calculated  by  means  of  the 
equation 


a(x) 


a 
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x da  dE 
0 dE  dx  dX  ’ 


(3.7) 


where  a is  the  K-shell  ionization  cross  section  for  the  incident  beam 

o 

energy.  The  average  depth  of  x-ray  emission  through  the  front  of  the 
target  is  then  given  by 


_ xa(x)T(x)dx 

x - . (3.8) 

/ 0 a(x)T(x)dx 

The  binary  encounter  model  results  of  Garcia**7  were  used  to  obtain  o 

o 

and  da/dE,  and  the  dE/dx  values  were  derived  from  the  stopping-power 
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tables  of  Northcliffe  and  Schilling.  Also,  it  was  necessary  to 
determine  the  average  projectile  energy  for  x-ray  emission.  The  average 
projectile  energy  for  the  emission  of  K x-rays  through  the  front  of  the 
target  was  calculated  from  the  relationship 


_ /!:  E(x)T(x)a(x)dx 

E --7 , (3.9) 

fQ  T(x)o(x)dx 


where  E(x)  - E - 4^  dx  . 

o 0 dx 


(3.10) 
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In  Eo.  (3.10),  Is  the  incident  beam  energy.  The  effective  target 

thicknesses  and  the  average  projectile  energies  were  calculated  using  a 
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computer  program,  ORION. 

For  the  purpose  of  making  detailed  comparisons  of  the  Ka  x-ray 

satellite  intensity  distributions  in  this  study,  it  is  convenient  to  use 

a parameter  which  shall  be  referred  to  as  the  apparent  average  L-vacancy 

fraction,  p This  parameter,  p , is  defined  as 

L L 
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where  n is  the  number  of  L-shell  vacancies,  N is  the  number  of  L-shell 

Li 

electrons  in  the  ground-state  atom,  and  f is  the  fraction  of  the  total 

n 

Ka  x-ray  yield  contained  in  the  nth  satellite  peak.  The  value,  p , does 

Li 

not  directly  represent  the  average  fraction  of  L-shell  vacancies  at  the 
time  of  Ka  x-ray  emission  since  the  fluorescence  yields  differ  for  the 
various  vacancy  configurations.  However,  this  parameter  has  been  chosen 
since  it  has  been  shown^’^  that  the  intensity  pattern  of  the  Ka 
satellite  lines  can  be  well  represented  by  binomial  distributions  in 
which  the  relative  intensity  of  each  satellite  peak  is  approximately 
given  by  the  formula 

I N.  N - n 

f - j2 — - (J ) P?  (1  - PL)  . (3.12) 

n Itot  n L L 
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where  ( ) ■ -rr rr — r—  , I is  the  number  of  x-rays  detected  in  the 

n (N  - n)ln!  * n J 

L 

nth  satellite  peak,  and  I t is  the  total  number  of  x-rays  detected. 
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CHAPTER  IV 
RESULTS 

Ka  x-ray  spectra  obtained  in  the  present  work  illustrating  the 
dependence  of  the  p^  values  on  the  projectile  atomic  number  are  shown 
in  Figs.  8 and  9.  As  the  atomic  number  of  the  projectile  is  increased, 
the  intensity  distribution  shifts  to  the  higher  order  satellites 
yielding  higher  p^  values.  This  effect  has  been  noted  previously  using 
thick^’^  and  thin"^  targets.  In  the  neon  on  silicon  monoxide  spectrum 
in  Fig.  8,  seven  peaks  are  discernible  whereas  in  the  aluminum  spectra 
in  Fig.  9 using  carbon,  oxygen,  and  neon  ions,  only  six  are  discernible. 
The  rapid  attenuation  of  the  sixth  and  seventh  satellite  peaks  in  the 
aluminum  spectra  is  attributed  to  the  K absorption  edge  occuring  at 
1559.6  eV. 

The  dependence  of  p on  the  velocity  of  the  projectile  is  shown  in 
L 

Figs.  10  and  11.  Fig.  10  illustrates  the  velocity  dependence  with  light 

(helium)  ions  and  Fig.  11  with  heavier  (oxygen)  ions.  In  both 

situations,  as  the  projectile  energy  is  increased,  the  intensity 

distribution  shifts  to  the  lower  order  satellites  or  to  a lesser  degree 

49 

of  L-shell  ionization.  A study  by  Li  jit  a^.  showed  that  for  deuterons 

and  helium  ions  on  calcium,  the  fraction  of  K(5  x-rays  emitted  in  the 
presence  of  L-shell  vacancies  reached  a maximum  at  a projectile-to- 
L-shell  electron  velocity  ratio  of  about  1.  The  results  of  Richard  et^ 
al. ^ studying  the  Ka  x-rays  from  helium  ions  on  aluminum  showed  that 
maximum  L-shell  ionization  in  the  presence  of  a single  K-shell  vacancy 


FIG.  9.  Aluminum  Ka  x-ray  spectra  produced  by  He,  C,  0,  and  Ne 
ions. 


I 


1 

FIG.  10.  Ka  x-ray  satellite  dependence  on  the  projectile  energy 
produced  by  light  (He)  ions. 


RELATIVE  COUNTS/CHANNEL 


35 


FIG.  11.  Ka  x-ray  satellite  dependence  on  the  projectile  energy 
produced  by  heavy  (0)  ions. 


1 
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occured  at  a velocity  ratio  of  approximately  1.25.  The  investigation  by 
54 

Knudson  et^  al.  of  Ka  x-rays  from  low  energy  neon  ions  on  aluminum 
indicated  that  maximum  L-shell  ionization  in  the  presence  of  one  K-shell 
vacancy  occurs  above  a velocity  ratio  of  2.  From  these  results,  it  is 
evident  that  the  velocity  ratio  for  maximum  L-shell  ionization  in  the 
presence  of  a K-shell  vacancy  is  dependent  upon  the  atomic  number  of  the 
projectile.  In  this  study,  ions  having  energies  above  the  velocity  ratio 
for  maximum  L-shell  ionization  were  used  and  thus,  in  Fig.  10,  the  p 

Li 

values  decrease  with  increasing  projectile  energy.  This  is  qualitative- 
ly the  behavior  expected  on  the  basis  of  the  velocity  dependence  of  the 
cross  section  for  L-shell  ionization. 

In  Figs.  10  and  11,  the  dependence  of  the  peak  to  background  ratio 
on  the  projectile  energy  may  also  be  observed.  The  background  rapidly 
increases  as  the  projectile  energy  increases.  In  Fig.  10,  a higher  rate 
of  increase  of  the  background  in  the  silicon  monoxide  spectra  than  in 
the  aluminum  spectra  is  seen.  The  more  rapid  rate  of  increase  is  a 
result  of  the  lower  Coulomb  barrier  for  helium  ions  on  the  oxygen  atoms 
in  the  silicon  monoxide  target  compared  with  the  Coulomb  barrier  of 
helium  ions  on  aluminum  atoms.  Nuclear  reactions  can  occur  at  lower 
energies  in  the  silicon  monoxide  target  and  thus,  the  more  rapid  rate  of 
increase  of  the  background  with  increasing  projectile  energy.  (A 
further  discussion  of  peak  to  background  measurements  can  be  found  in 
Appendix  B. ) 

The  measured  Ka  x-ray  relative  intensities  are  listed  in  Table  I 
for  helium,  carbon,  and  neon  ions  on  aluminum  and  silicon.  Also  listed 
in  this  table  are  the  incident  projectile  energies,  the  average 


39 


TABLE  I.  Relative  Ka  x-ray  intenaitlea  and  values. 


Projectile 

and 

t 

n 

Compound 

n*0 

n-1 

ri*2 

n-3 

n-4 

n-5 

n-6 

NOM 

E 

PL 

5.42  MeV  He 


A1 

0.580 

0.336 

0.084 

5 

4.86 

0.06310.001 

S10 

0.550 

0.371 

0.079 

6 

5.04 

0.066±0.001 

6.64  MeV  He 

A1 

0.643 

0.285 

0.064 

3 

6.14 

0.05210.002 

S10 

0.630 

0.317 

0.053 

3 

6.61 

0.053*0.001 

10.40  MeV  He 

A1 

0.761 

0.208 

0.032 

3 

10.03 

0.034+0.002 

S10 

0.753 

0.219 

0.029 

3 

10.05 

0.03510.003 

23.20  MeV  He 

A1 

0.865 

0.135 

5 

22.94 

0.017+0.0005 

SiO 

0.849 

0.151 

5 

23.03 

0.019+0.0007 

40.60  MeV  He 

A3. 

0.895 

0.105 

5 

40.45 

0.013+0.0005 

SiO 

0.896 

0.104 

5 

40.48 

0.013+0.001 

12.00  MeV  C 

A1 

0.029 

0.120 

0.337 

0.292 

0.161 

0.061 

3 

9.35 

0.32710.002 

SiO 

0.027 

0.151 

0.308 

0.304 

0.138 

0.073 

3 

9.81 

0.32410.0002 

22.00  MeV  C 

A1 

0.058 

0.181 

0.377 

0.251 

0.105 

0.029 

3 

18.51 

0.282i0.001 

SiO 

0.061 

0.238 

0.377 

0.235 

0.087 

3 

18.77 

0.25610.003 

26.90  MeV  C 

A1 

0.078 

0.219 

0.396 

0.223 

0.084 

3 

23.58 

0. 252lO.  001 

SiO 

0.076 

0.274 

0.352 

0.277 

0.071 

3 

24.83 

0.24310.001 

42.9  MeV  C 

A1 

0.177 

0.318 

0.344 

0.129 

0.033 

4 

40.31 

0.190t0.G02 

SIO 

0.172 

0.373 

0.304 

0.127 

0.024 

7 

41.23 

0.18210.001 

74.1  MeV  C 

Ai 

0.316 

0.360 

0.324 

4 

72.27 

0.12610.003 

SiO 

0.350 

0.438 

0.213 

4 

72.88 

0.10810.0006 

23.4  MeV  Me 

Al 

0.032 

0.081 

0.254 

0.289 

0.220 

0.124 

5 

17.81 

0.37010.002 

SiO 

0.021 

0.092 

0.221 

0.281 

0.189 

0.143 

0.053 

3 

19.80 

0.39510.001 

46.1  MeV  He 

Al 

0.055 

0.089 

0.278 

0.279 

0.196 

0.103 

3 

34.05 

0.34810.0007 

SiO 

0.040 

0.106 

0.247 

0.286 

0.172 

0.112 

0.037 

3 

36.56 

0.36610.0002 

117.0  MeV  Ne 

Al 

0.142 

0.226 

0.332 

0.216 

0.070 

0.011 

3 

112.62 

0.23510.002 

SiO 

0.121 

0.253 

0.330 

0.211 

0.069 

0.017 

3 

114 . 11 

0.238i0.001 

22.00  MeV  C 

Si (Thick) 

0.057 

0.225 

0.347 

0.270 

0.102 

3 

17.29 

0.26710.001 

SiO(Thick) 

0.058 

0.230 

0.357 

0.260 

0.095 

3 

17.81 

0.26310.001 

S102 (Thick) 

0.062 

0.234 

0.366 

0.253 

0.084 

3 

18.04 

0.25810.001 
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projectile  energy  for  x-ray  emission  and  detection  (E) , the  number  of 

measurements  taken  for  each  compound  or  element  at  the  given  projectile 

energy  (NOM) , and  the  p values.  The  f values  (the  ratio  of  the  nth 

L n 

satellite  peak  intensity  to  the  sum  of  the  intensities  of  all  the  Ka 
peaks)  have  been  corrected  for  absorption  in  the  target  and  detector 
window,  and  for  detection  efficiency.  The  p^  values  were  determined 
directly  from  the  f^  values  and  the  indicated  errors  are  standard 
deviations.  Table  I shows  numerically  the  shift  in  the  intensity 
distribution  to  the  lower  order  satellites  with  increasing  projectile 
energy. 

Also  listed  in  Table  I are  the  relative  intensities  and  p values 

L« 

for  22.0  MeV  (incident  energy)  carbon  ions  on  thick  silicon,  silicon 
monoxide,  and  silicon  dioxide  targets.  The  variations  of  the  intensity 
distributions  and  the  p values  from  one  target  to  the  next  indicate 

Li 

differing  average  L-shell  vacancies  at  the  time  of  Ka  x-ray  emission. 

The  differing  average  L-shell  vacancies  at  the  time  of  Ka  x-ray  emission 

may  be  associated  with  the  vacancy  production  process  which  occurs 

during  collision,  or  with  the  deexcitation  process.  Earlier  investiga- 
39  47 

tions  ’ have  indicated  that  at  the  beam  energies  used  in  these 

measurements,  it  is  quite  improbable  that  the  chemical  environment  of  the 

target  atom  could  influence  the  vacancy  distribution  produced  during  the 

collisions.  The  studies  have  also  indicated  that  these  variations  may 

be  attributed  to  intra-  and  interatomic  transitions  occuring  prior  to  Ka 

x-ray  emission.  These  transitions  can  contribute  in  a major  way  to  K- 

and  L-shell  vacancy  filling  in  multiply  ionized  atoms. 

A plot  of  the  p values  as  a function  of  the  velocity  ratio  for  the 
Li 


target-projectile  combinations  listed  in  Table  I is  shown  in  Fig.  12. 
The  velocity  ratio  was  calculated  from  the  expression 


vi  m E.  , 

_1  . (-»-!)*, 

VL  M UL 


(4.1) 


where  is  the  average  projectile  velocity  for  x-ray  emission  with 

energy  E , V is  the  weighted  average  of  the  velocities  of  the  L-shell 
X Li 

electrons,  m^  is  the  electron  rest  mass,  M is  the  mass  of  the  projectile, 
and  U is  the  weighted  average  of  the  L-shell  electron  binding  energies. 

Li 

The  L-shell  binding  energies  used  in  this  calculation  were  taken  from  a 

72 

tabulation  by  Bearden  and  Burr. 

The  dependence  of  p on  the  energy  and  atomic  number  of  the 

Lj 

projectile  and  on  the  atomic  number  of  the  target  atom  is  shown  in 
Fig.  12.  As  the  projectile  energy  increases,  the  p values  decrease. 

L 

As  the  projectile  atomic  number  increases,  the  p values  increase. 

L 

Finally,  as  the  atomic  number  of  the  target  increases,  the  p values 

Li 

decrease.  Also  shown  in  Fig.  12  is  the  flattening  of  the  curves  with 

increasing  projectile  atomic  number  and  the  apparent  convergence  of  the 

heavy  ion  data  at  low  energies  as  the  region  of  maximum  p is  approached. 

Li 

This  is  suggestive  of  a saturation  effect  in  which  the  x-ray  spectra  and 

subsequent  p values  reach  a point  where  they  no  longer  reflect  the 
L 

initial  vacancy  distribution  as  a result  of  L-shell  vacancy  filling. 

The  p values  for  helium  ions  on  aluminum  corrected  for  fluorescence 
Li 

yield  are  shown  in  Fig.  13.  The  low  energy  data,  indicated  by  open 
circles,  and  the  fluorescence  yields  were  obtained  from  a previous  study 


I 
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FIG.  12.  Variation  of  with  projectile  energy  for  A1  and  SiO 
targets. 


Fluorescence  Yield  Corrections 
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by  Richard  £t  al.  The  fluorescence  corrections  were  made  according 
to  the  expression 


I' 

n 


I 

n 


0) 

KL 


n 


(4.2) 


where  I is  the  intensity  of  the  nth  satellite,  I'  is  the  corrected 
n n 

intensity  of  the  nth  satellite,  and  u)  is  the  fluorescence  yield  for 

KLn 

the  vacancy  configuration  with  one  K-shell  and  n L-shell  vacancies. 

The  fluorescence  yields  used  for  the  various  configurations  were  0.041, 
0.042,  and  0.044  corresponding  to  n=0,  1,  and  2 respectively."*"*  An 
assumption  made  by  Richard  et_  al^.  "*"*  in  the  determination  of  the 
fluorescence  yield  was  that  the  L-shell  vacancies  were  in  the  and 
Ljjj.  subshells  (2p  electrons).  This  assumption  is  not  entirely  valid 
and  some  errors  will  be  introduced  into  the  fluorescence  yield 
corrections.  It  is  apparent  that  the  fluorescence  yield  corrections  do 
not  change  the  features  of  the  curve  in  a significant  way  for  the  low 
states  of  ionization  involved  for  this  projectile-target  combination 
(an  average  of  1.2  L-shell  vacancies  in  the  atoms  of  the  highest  state 
of  ionization  shown  in  Fig.  13) . 


CHAPTER.  V 


COMPARISON  WITH  THEORY 
A.  General  Formulations 


The  collision  of  an  energetic  projectile  ion  with  a stationary 
target  atom  generally  results  in  the  transfer  of  some  of  the  energy  of 
the  projectile  to  the  target  atom.  Frequently,  tills  energy  transfer 
results  in  the  ionization  of  one  or  more  inner-shell  electrons.  Most 
simple  models  of  inner-shell  ionization  begin  with  the  assumption  that 
the  production  of  an  inner-shell  vacancy  occurs  as  a result  of  the 
direct  Coulomb  interaction  of  the  incident  particle  with  the  bound 
electron.  Three  such  models  - the  plane  wave  Born  approximation 
(FVBA) ^3,74  impulse  or  binary--  encounter  approximation 

(BF.A) , ^ ^ and  the  semi-classical  approximation  (SCA)7^  - will 

be  briefly  discussed  here  and  the  results  of  two,  the  BE A and  SCA,  will 
be  compared  with  experimental  data. 

The  PWBA.  and  BEA  are  high  energy  formulations  and  are  expected  to 

be  valid  for  incident  particles  with  energies  much  greater  than  the 

binding  energy  cf  the  atomi a electron.  The  high  energy  criterion  for 

74  79 

these  approximations  is  generally  given  as  ' 


7lV' 

“ hv. 


« 1 


(5.1) 


where  is  the  atomic  number  of  the  p>ojccti.le,  Is  the  atomic  nnmbe 
of  the  target,  e is  the  charge  of  the  electron,  h is  tin  Planck  constan 
divided  by  2n,  and  v^  is  l he  projectile  velocity. 
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The  first  order  plane  wave  Born  approximation  assumes  that;  1)  the 

projectile  acts  as  a point  charge  (the  effect  of  the  atomic  structure 

of  the  projectile  on  vacancy  production  is  negligible),  2)  plane  waves 

can  be  used  to  describe  the  incident  and  scattered  particle  over  all 

space,  and  3)  the  electron  transition  is  from  the  bound  state  to  the 

continuum  with  the  other  electrons  remaining  in  their  initial  states. 

The  ionization  cross  sections  for  helium  on  aluminum  predicted  by  the 

74 

PWBA  can  be  calculated  by 
8irn  z2a2 

°±m— — -vw  • (5*2) 

Zeffni 

where  n^  is  the  number  of  electrons  in  the  i-th  shell,  z^  is  the 

projectile  atomic  number,  a is  the  Bohr  radius,  z is  the  charge  of 

o ef  f 

the  target  nucleus  corrected  for  screening  by  the  inner-shell  electrons, 

and  f i are  form  factors  calculated  using  non- relativistic  hydrogenic 

wave  functions.  These  form  factors  have  heen  calculated  and  tabulated 

80  81 

as  functions  of  0^  and  by  Khandelwal  et^  al.  and  Choi  ^ al.  for 
the  K-  and  L-shells,  and  were  used  to  determine  the  PWBA  ionization 
cross  sections  shown  in  Figs.  14  and  15.  The  quantity,  6^,  the  ratio 
of  the  true  binding  energy  of  the  electron  to  that  predicted  by  a 


screened  hydrogenic  wave  function,  is  determined  by  the  expression 


I 
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FIG.  14.  Theoretical  K-shell  ionization  cross  sections. 


FIG.  15.  Theoretical  L-shell  ionization  cross  sections. 
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quantity,  ru,  is  given  by  the  expression 


'ef f e 


hv,  m E, 

^ 2 ' M.  2 


1 ZeffR 


(5.4) 


where  m is  the  rest  mass  of  the  electron  and  M,  and  E„  are  the  mass  and 
o 11 

energy  of  the  projectile  respectively. 

In  the  impulse  or  binary-encounter  approximation,  ionization  is 

treated  classically  as  a collision  between  the  projectile  and  a free 

electron.  The  result  of  the  collision  is  a direct  energy  transfer  or 

exci  ange  between  the  charged  particle  and  the  electron.  All  other 

interactions  are  neglected  and  the  target  nucleus  exists  simply  to 

8 2 

establish  a momentum  distribution  for  the  electron.  Thomas  and 
83 

Garcia  have  adjusted  this  model  to  take  into  account  the  effects  of 

nuclear  repulsion  between  the  projectile  and  the  target  nucleus.  Garcia 
67  83 

et  al.  * have  shown  that  if  hydrogenic  velocity  distributions  are 
used  for  the  exchange  of  energy  between  the  two  particles,  the  BEA 
ionization  cross  sections  for  the  i-th  shell  obeys  a scaling  law 


2 2 R1 
Ui  = Zlf(^> 


for 


M 


<<  1. 


(5.5) 


Thus,  the  product  of  the  binding  energy  squared  times  the  ionization 

2 

cross  section  divided  by  z^  is  a universal  function  of  the  incident 
projectile  energy  expressed  in  units  of  the  binding  energy.  Using  a 
semi-relativistic  velocity  ratio. 
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in  place  of  the  classical  velocity  ratio,  and  performing  a linear 

83 

interpolation  in  Garcia's  tables,  the  BEA  cross  sections  for  helium 
ions  on  aluminum  were  determined  for  the  K-shell  (Fig.  14).  The  BEA 
L-shell  ionization  cross  sections  for  helium  ions  on  aluminum  were 
obtained  from  a tabulation  of  similar  calculations  by  Hansen"^  and  are 
shown  in  Fig.  15. 

The  impact  parameter  or  semi-classical  approximation  was  formulated 
78 

by  Bang  and  Hansteen  to  treat  ion-atom  collisions  that  did  not  satisfy 
the  high  energy  criterion  of  Eq.  (5.1).  This  treatment  is  identical  to 
the  PWBA  with  the  exception  that  corrections  are  made  for  the  deflection 
of  the  projectile  near  the  target  nucleus  as  a result  of  Coulomb 
repulsion.  The  repulsion  prevents  the  projectile  from  penetrating  too 
closely  to  the  target  nucleus  and  thus,  cross  sections  lower  than  the 
PWBA  predictions  can  be  expected.  The  SCA  predictions  of  the  K-  and 
L-shell  ionization  cross  sections  for  helium  ions  on  aluminum  are 
given  in  Figs.  14  and  15  and  were  determined  by^ 
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where  g is  equal  to  unity  for  the  K-shell  and  Lj-subshell,  2/3  for  the 
Ljj-subshell,  and  4/3  for  the  Lj^-subshell,  and  Zgff  Is  the  screened 
nuclear  charge  given  by 


'«££  ‘ z2  ' S ’ 


(5.8) 


where  S is  the  Slater  screening  coefficient  for  the  shell  being  ionized. 

The  term,  0^,  is  identical  to  Eq.  (5.3).  The  values  of  F^  were  obtained 

84 

from  a tabulation  by  Kocbach  as  a function  of  x where  x is  determined 
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by 


226i 

“ir* 


(5.9) 


The  shapes  and  magnitudes  of  the  PWBA,  BEA,  and  SCA  predictions 
for  K-shell  ionization  cross  sections  shown  in  Fig.  14  are  comparable. 

An  unexpected  feature  in  Fig.  14  is  the  difference  between  the  PWBA 
and  SCA  curves.  One  would  expect  the  SCA  predictions  to  be  lower  than 
the  PWBA  predictions  at  lower  energies,  and  comparable  at  higher 
energies.  In  Fig.  14,  the  SCA  predictions  are  slightly  higher  at 
lower  energies  and  considerably  lower  at  the  higher  energies. 

The  PWBA,  BEA,  and  SCA  differ  greatly  in  predicting  the  shape 
and  magnitude  of  the  dependence  of  L-shell  ionization  on  the  velocity 
ratio.  In  Fig.  15,  the  PWBA  and  BEA  predictions  of  the  shape  of  the 
dependence  of  the  ionization  cross  section  on  the  velocity  ratio  differ 
considerably,  however,  the  overall  trend  of  the  curves  are  similar. 

The  SCA  predictions,  on  the  other  hand,  are  as  much  as  an  order  of 
magnitude  lower  than  the  PWBA  and  BEA  predictions  and  show  no  indication 
of  a maximum  at  the  lowest  velocities  for  which  tabulated  values  of  F^ 
are  avaibable. 

B.  BEA  L-Shell  Ionization  Probabilities 

The  BEA  probability  for  ionization  of  an  L-shell  electron  as  a 
function  of  impact  parameter  was  obtained  using  a computer  program, 
BEACON,  by  Hansen. ^ In  the  Hansen  impact  parameter  representation  of 
the  BEA,  the  ionization  cross  section  of  the  i-th  shell,  a^,  is  related 
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to  the  total  ionization  probability,  P^,  by 

°i  “ V Pi(b)db  * (5.10) 

The  probability  of  ionization  per  electron,  p^,  is  related  to  the  total 
probability  by 

Pi(b)  = 1 - [1-Pl(b)]n  , (5.11) 

where  n is  the  number  of  electrons  in  the  i-th  shell. 

At  selected  energies,  the  probabilities  of  ionization  per  electron 
for  helium  ions  on  aluminum  were  obtained  as  a function  of  impact 
parameter  for  the  K-  and  L-shells  (2s  and  2p  electrons) . The  L-shell 
probabilities  were  summed  at  given  impact  parameters  and  divided  by  the 
number  of  L-shell  electrons  to  obtain  the  average  probability  per 
electron  as  a function  of  impact  parameter. 

C.  SCA  L-Shell  Ionization  Probabilities 

The  SCA  probability  for  ionization  of  an  L-shell  electron  as  a 
function  of  impact  parameter  for  a given  shell  or  subshell,  i,  and 
characterized  by  the  hydrogenic  quantum  numbers,  were  determined 

by 


2j.+l  1 

p(b)  “ Vn  # (04»XJ — I . (X,,B.)  , 

i’*i  1 1 21  +1  z^S  i,£,i  1 1 


(5.12) 


where  j is  the  angular  momentum  of  the  shell  or  subshell,  0^  is 
described  in  Eq.  (5.3),  and  z^  is  described  in  Eq.  (5.8).  The  values 
of  were  calculated  by 
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(5.13) 


where  here  is  the  principle  quantum  number  of  shell,  i,  and  is 
the  projectile  energy  in  MeV.  The  values  of  3^  were  determined  by  the 
expression 
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The  values  for  p (0.»X.,)  and  I „ (y. >6.)  were  then  obtained  from 

n^,£^  i i n^,£.^  1 1 
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a tabulation  by  Hansteen  eit  al.  as  functions  of  0^  and  x^»  and  and 

6^  respectively  and  substituted  into  Eq.  (5.12)  to  obtain  the  ionization 

probabilities  as  a function  of  impact  parameter.  Using  this  procedure, 

the  ionization  probabilities  for  protons  on  aluminum  were  determined. 

2 

Simple  z scaling  was  used  to  predict  the  ionization  probabilities  for 
helium  ions  on  aluminum.  At  each  energy,  the  probabilities  as  a 
function  of  impact  parameter  were  obtained  for  the  K-shell  and  the  L^-, 
L^-,  and  L^. ^-subshells.  The  L-subshell  probabilities  were  summed  at 
the  given  impact  parameter  and  divided  by  the  number  of  L-shell 
electrons  to  obtain  the  average  probability  per  electron  as  a function 
of  impact  parameter. 


D.  Comparison  of  Theoretical  Predictions 
with  Experimental  Data 

Polynomial  fits  using  a computer  program,  POLLY,  were  made  to  the 
calculated  points  of  the  probabilities  of  ionization  per  electron  as  a 
function  of  impact  parameter  for  each  shell  (K-  and  L-shells)  to  obtain 


V 

I 


56 


the  functions,  PR(b)  and  P^Cb).  The  functions  PR(b)  and  ?L(b)  are  then 
the  probability  per  electron  of  ionizing  one  K-shell  electron  and  the 
average  (2s  and  2p)  probability  per  electron  of  ionizing  one  L-shell 
electron  respectively.  Using  P (b)  and  P (b) , the  theoretical 

K.  Li 

ionization  cross  sections  for  producing  a state  having  one  K-shell 
and  n L-shell  vacancies  were  determined  by 

o n = 27T/y(J)P  (b)(l-P  (b))(®)P  (b)n(l-P  (b))8“nbdb.  (5.15) 

n OIK  K nL  L 

IUj 

The  integration  of  Eq.  (5.15)  over  impact  parameters  was  accomplished 

by  the  Gaussian  numerical  integration  method  described  in  Appendix  A. 

The  upper  limit  of  the  integral  was  determined  by  the  K-shell 

probabilities  and  was  selected  such  that  the  K-shell  probability  at  the 

4 

upper  limit  was  more  than  a factor  of  10  smaller  than  the  maximum 
K-shell  probability  at  that  energy.  The  lower  limit  of  the  integral 
was  selected  to  be  zero. 

The  p value  for  each  energy  was  then  calculated  from  the 
Li 

ionization  cross  sections  by 


i!  KLn 

pL  " 8 E,” 

n“l  tot 


(5.16) 


where  o is  the  sum  of  the  ionization  cross  sections  of  all  the 
tot 

vacancy  configurations. 

The  results  of  the  BEA  and  SCA  predictions  of  pL  on  the  projectile- 

to-L-shell  electron  velocity  ratio  are  shown  in  Fig.  16.  Also  shown  in 

Fig.  16  are  the  experimental  p values  from  helium  and  neon  ions 

L 

incident  on  aluminum.  In  Fig.  16(a),  the  low  energy  neon  on  aluminum 
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data  were  obtained  from  a study  by  Knudson  et^  al.  and  are  indicated 
by  open  circles.  The  low  energy  helium  on  aluminum  data  in  Fig.  16(b) 
indicated  by  open  circles  were  obtained  from  a study  by  Richard  et  al.'* 
The  neon  p values  show  a very  broad  maximum  centered  about  a velocity 

L 

ratio  of  2.3  and  decrease  very  slowly  with  increasing  velocity  ratio. 

The  experimental  helium  p values  in  Fig.  16(b),  on  the  other  hand, 

L 

reach  a very  sharp  maximum  at  a velocity  ratio  of  about  1.25  and 
decrease  much  more  rapidly  than  the  neon  p values  with  increasing 

L< 

velocity  ratio.  The  data  for  helium  on  aluminum  converge  at  high 

velocities  with  the  shakeoff  limit  observed  in  photoionization  and 

85 

electron  bombardment  studies  (p^=0.013).  Very  good  agreement  is 
exhibited  here  between  the  two  sets  of  measurements  for  helium  and  neon 
ions  on  aluminum. 

The  BEA  predictions  for  helium  on  aluminum  are  a factor  of  5 to  10 
times  larger  in  absolute  magnitude  than  the  experimental  values. 
Although  the  shape  of  the  velocity  dependence  is  not  well  described  by 
the  BEA,  a maximum  is  predicted  around  a velocity  ratio  of  2.  The  BEA 
results  for  helium  on  aluminum  are  surprisingly  similar  to  the  experi- 
mental results  for  neon  on  aluminum.  The  SCA  predictions  more  closely 
represent  the  slope  of  the  experimental  p values  in  the  high  velocity 

Lj 

region,  however,  they  are  a factor  of  about  3 too  low.  In  the  low 
velocity  region,  greater  divergence  from  the  experimental  data  is  seen. 
The  SCA  fails  to  predict  the  maximum  in  the  region  where  this  feature 
is  found  experimentally. 

The  helium  on  aluminum  combination  was  used  to  test  the  applica- 
bility of  the  BEA  and  SCA  to  the  description  of  multiple  ionization 
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because  of  the  low  states  of  ionization  resulting  from  these  impacts. 
These  low  states  of  ionization  (an  average  of  1.2  L-shell  vacancies  for 
the  highest  state  of  ionization  in  Fig.  16(b))  are  not  expected  to  be 
greatly  influenced  by  variations  in  the  fluorescence  yields  for 
different  vacancy  configurations  (as  shown  in  Fig.  13)  or  by  fast 
rearrangement . 

Similar  calculations  were  made  for  the  heavier  projectiles  but 
these  proved  to  be  meaningless.  The  BEA  and  SCA  greatly  over  predicted 
L-shell  ionization  to  the  extent  that  the  probabilities  of  ionizing  one 
L-shell  electron  was  essentially  equal  -o  unity  at  all  energies. 
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CHAPTER  VI 
CONCLUSIONS 

The  results  of  this  study  have  provided  a comprehensive  picture  of 
the  dependence  of  the  aluminum  and  silicon  Ka  x-ray  satellite  structure 
on  the  atomic  number  and  energy  of  the  projectile  at  high  energies.  The 
energies  used  in  this  study  were  above  the  energy  region  for  maximum 
L-shell  ionization  and  consequently,  decreasing  p values  were  observed 

L 

with  increasing  projectile  energy.  The  dependence  was  also  observed  to 
flatten  with  increasing  projectile  atomic  number  as  the  region  of  the 
PL  maximum  was  approached.  This  flattening  is  suggestive  of  a saturation 
effect  in  which  the  x-ray  satellite  spectra  of  heavy  ions  no  longer 
reflect  the  initial  vacancy  distribution.  The  initial  vacancy  distri- 
bution may  be  altered  by  L-shell  vacancy  filling  in  the  highly  ionized 
atom  prior  to  x-ray  emission  by  intra-  and  interatomic  processes.  The 
silicon,  silicon  monoxide,  and  silicon  dioxide  results  indicate  that 
fast  rearrangement  does  occur  to  some  extent. 

The  data  from  this  study  were  found  to  be  in  excellent  agreement 

with  the  results  of  previous  studies  by  Richard  £t  al. ^ and  Knudson  jet 

54 

al.  Together,  they  provide  a very  good  picture  of  the  dependence  of 
the  p^  values  on  the  energy  of  the  projectile  for  helium  on  aluminum 

and  neon  on  aluminum  over  a broad  energy  region. 

The  comparison  of  the  experimental  data  with  theoretical  predic- 
tions were  rather  discouraging.  The  BEA  greatly  overpredicted  the 

magnitude  of  the  p values  and  failed  to  accurately  predict  the  shape 
Li 
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of  the  dependence.  The  SCA  predictions  more  closely  represented  the 
slope  of  the  experimental  p values  at  the  higher  velocities  but  were 

L< 

too  low  in  magnitude.  At  the  lower  energies,  greater  divergence  from 
the  experimental  values  were  observed.  The  results  of  the  present 
investigation  clearly  demonstrate  that  a much  more  sophisticated 
theoretical  treatment  of  multiple  inner-shell  ionization  will  be 
required  for  a full  understanding  of  the  velocity  dependence  of  this 
process. 
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APPENDIX  A 
GAUSSIAN  INTEGRATION 


The  integral,  / f(x)dx,  can  be  computed  by  Gauss'  formula  often 
expressed  by  the  equations 


f (z)dz 

(A.  1) 

n 

I1wif(z1) 

(A. 2) 

n/2 

= w f(0)  + I w [f(z  ) + f(-z  )]  , 
i=l  1 1 


(A. 3) 


where  n is  the  number  of  interpolation  points,  are  the  weights,  and 
z^  are  the  abscissas.  The  interval,  [A,B],  can  be  transformed  to 


Let  x = cz 

+ d, 

then. 

dx  = cdz 

(A. 4) 

A = -c  + d 

when 

z=-l , and 

(A. 5) 

B = c + d 

when 

N 

II 

h-» 

(A. 6) 

Combining  Eqs.  (A. 5)  and  (A. 6),  the  following  relations  are  obtained: 


c 


From  Eq. 


B - A . , A + B 

— r — and  d * — - — 


(A. 4),  if  dx  = cdz,  then 
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g(x)  = cg(cz  + d)  and 

f ( z ) = eg  (cz  + d).  (A. 7) 

Equations  (A.l)  to  (A. 3)  then  become 

B m 

/ g(x)dx  = E cw  g (cz  + d),  (A. 8) 

A i=l  1 

where  in  is  a factor  that  determines  the  number  of  interpolation  points, 
n. 

The  n-point  Gaussian  integration  routine  listed  at  the  end  of  this 
appendix  was  written  to  calculate  the  integral, 

Q n = 2^A(l)PK(b)(1"PK(b)(n)(?L(b))n(1_\(b))8  nbdb  ’ (A’9) 

KL 

as  described  in  Chapter  V.  P (b)  and  P (b)  can  be  any  polynomial 

K Li 

function  up  to  and  including  a fifth  order  polynomial.  The  number  of 
interpolation  points  can  be  selected.  The  data  card  formats  are  as 
follows : 


Card  1 - title  card 


Columns 

Format 

Variable 

Identity 

3-80 

A6 

Any  title  or  description 

Card  2 - 

coefficients  for  the 

5 4 3 

equation,  S6x  + s5x  + S4x 

+ S3x2  + 

S2x  + SI 

Columns 

Format 

Variable 

Identity 

1-10 

F10.6 

SI 

Coefficients 

11-20 

F10.6 

S2 

for 

21-30 

F10.6 

S3 

pK(h) 
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I 

t 


Columns 

Format 

Variable 

Identity 

31-40 

F10.6 

S4 

41-50 

F10.6 

S5 

51-60 

F10.6 

S6 

Card  3 - 

5 4 3 

coefficients  for  the  equation,  T6x  + T5x  + T4x 

T3x2  + T2x 

+ T1 

Columns 

Format 

Variable 

Identity 

1-10 

F10.6 

T1 

Coefficients 

11-20 

F10.6 

T2 

for 

21-30 

F10.6 

T3 

31-40 

F10.6 

T4 

Vb> 

41-50 

F10.6 

T5 

51-60 

F10.6 

T6 

Card  4 - 

integral  limits  and 

interpolation  points 

Columns 

Format 

Variable 

Identity 

1-10 

F10.6 

A 

Upper  limit  of  integral 

11-20 

F10.6 

B 

Lower  limit  of  integral 

21-22 

12 

m 

Factor  to  determine  number 
of  interpolation  points 

The  possible  values  for  m corresponding  to  the  number  of  interpolation 
points  are  as  follows: 

Values  of  m 2 3 4 5 6 10  15  20  40 

Number  of  interpolation  points,  nl2233  5 8 10  20 

The  ionization  cross  sections  in  the  output  are  not  absolute  cross 
sections.  The  probabilities  for  ionization  used  in  determining  PK(b) 


and  P (b)  were  normalized  for  use  in  the  polynomial  fitting  program 
L 
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i 

) 

I 

f 

and  the  normalized  functions  were  subsequently  used  in  the  ionization 
cross  section  calculations. 
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I JOB  76133C  BLAKE  SONOBE 

MBBOX  CYCLOTRON 

I B JOB  VERSION  5 HAS  CONTROL.  (IBSYS  VERSICN  13  3 
1 1 6 JOB 
MPFTC  MAIN 

DIMENSION  TITLEI13) 

COMMON  SI ,S2,S3,S4,S5,S6 
COMMON  T1 , T2, T3, T4,T5,T6 

1 IF  <EOF(5))  6C0.2.600 

2 REAO  15,8001  NC 
READ  (5,1001  TITLE 
WRITE  (6,200)  TITLE 
IF  (NC.E0.2)  CO  TO  3 
IF  (NC.EO.l)  GO  TO  4 
WRITE  (6,205)  NC 
GO  TO  600 

3 REAO  (5,103)  SI , S2 , S3, S4, S3 , S6  / 

WRITE  (6,202)  S6,S5,S4 

WRITE  (6,204)  S3 , S2 , SI 
A REAO  (5,103)  T l , T2 , T 3 , T4 , T5 , T6 
WRITE  (6,203)  T6,T5,T4 
WRITE  (6,204)  T3,T2,Tl 
REAO  (5,101)  A,B,M 
WRITE  (6,201)  A , B , M 

CALL  GAUSS  (A,8,M,SIG* NC ) / 

GO  TO  I 
C 

C ••••••INPUT  FORMATS. 

C 

100  FORMAT  1 2X , 13 A6 ) 

101  FORMAT  ( 2F 10*  6 , I 2 ) 

103  FORMAT  (6F10.6) 

800  FORMAT  (ID 

C 

C OUTPUT  FORMATS 

C 

200  FORMAT  ( IHI,10X,13A6) 

201  FORMAT  (1H-,10X,35HTHE  INPUT  VALUES  FOR  A, B , ANO  H ARE, 3X, F10.6, 

1 3X,F10.6,3X,3HAND,3X,12) 

202  FORMAT  ( 1H-,10X,29HTHE  INPUT  FUNCTION  FOR  FX1  I S , 3X , F 10 . 6 , 5H*X**5 , 
l 3H  ♦ ,F 10.6,5H*X**4,3M  ♦ , F 10 . 6 , 5H*X**3 , 2H  ♦ ) 

203  FORMAT  ( IH- , l OX , 29HTHE  INPUT  FUNCTION  FOR  FX2  IS, 3X.F10.6, 5H*X**5, 
1 3H  ♦ ,F10.6,5H*X**4,3M  ♦ , F 1 0. 6 , 5H*X • • 3 , 2H  ♦» 

204  FORMAT  ( 1H0,15X,F10.6,5H*X**2,JH  ♦ , F 10 . 6 , 2H#X , 3H  ♦ ,F10.6) 

205  FORMAT  I l HO , 1 OX , 4HNC  - , 1 1 , 2X , 1 9h 1 S AN  INVAL10  INPUT) 

600  STOP 

ENO 


GAUSSIAN  INTEGRATION 
INSTITUTE 
07/20/76 


tIBFTC  AA 

subroutine  gauss  ( A,B,MvS!G,NC) 

DIMENSION  NP0INTI9),  KEY(IO),  2154),  WEIGHTI54) 
OAT  A NPUI NT /2, 3, 4, 5, 6, 10, 15,20,40/ 

OAT  A KEY/l ,2 ,4,6,9,12,17 ,25,35,55/ 

OATA  NX/9/ 


DMA  l 

/O. 57735027 

,0.0  • 

,0.77659667 

2 COOl 

l 

0. 33990104 

,0.86113631 

,0.0 

,0.53866931 

ZC002 

2 

0.90617985 

,0.23861919 

,0.66120939 

,0.93266951 

ZC003 

3 

0.  14887434 

,0.63339539 

,0.67960957 

,0.86506337 

ZC006 

4 

0.97390653 

,0.0 

,0.20119609 

,0.39615135 

ZC003 

5 

0.57097217 

,0.72661773 

,0.86020668 

,0.93727339 

ZC006 

6 • 

0.90799252 

,0.07652652 

,0.22778585, 

1 0. 3 7 170609  , 

Z 0007 

70.51086  700,0.6  5604  J6B,0.  I 46  ill'll  .0. 839  U61I,0.'H  i 2 34  <.1,0. 96397 19  3,20008 

80. 99312860.0.  03877^42, 0.11608808, 0.19267758, 0.26815219, 0.34100407, /C009 

90.41377920.0. 48307580.0.54946713.0.61255389.0.67195668.0. 72731826.20010 
AO. 7 78 30565, 0.8246 122 3,0. 86575750,0. 90209  881, 0.9  3281281, 0.95771682, 2001 1 


00.97725995, 
DATA  HEIGHT 

0.99072626,0 

.99023771/ 

/ 1.0 

,0.88888889 

,0.55555556 

f 

10012 

WOOOl 

1 

0.65216516 

,0.36785686 

,0.56880889 

,0.67862867 

• 

WC002 

2 

0.23692688 

,0.66791396 

,0.36076157 

,0.17132669 

• 

WC003 

3 

0.29552622 

,0.26926672 

,0.21908636 

,0.16965135 

t 

W0004 

6 

0.C6667136 

,0.20267826 

,0.19863168 

,0.18616100 

• 

W0005 

5 

0.16626921 

,0.13957068 

,0.10715922 

,0.07036605 

• 

WC006 

6 

0.03075326 

,0.15275339 

,0.16917299, 

0. 16209611, 

t 

W0007 

70.1 3168866, 

0. 11819653,0 

.10193012,0 

.08327676,0. 

06267205,0.06060163 

. W0008 

80.01761401 .0.  07750595.0. 07703982. 0.0 76 110 36. 0.074 72 31 7.0. 07288658. 80009 
90. 0 706 1165, 0.06 791 205, 0.0648040 1,0.06130624, 0.05 74 397 7, 0.053 2 2 785, WOO  10 
AO. 04 869531 ,0. C4 38 709 1 , 0.0387 82 1 7 , 0. 03346020 , 0 .02793701 , 0. 022 24585, W001 1 

80.01642106.0. 01049828.0.00452128/  WC012 

IF  (NC.EQ.11  GO  TO  60 

IF  (NC.EQ.2)  GO  TO  50 
C 

C I ON  I 2 A T 1 ON  CALCULATION  OF  l L-SHELL  VACANCIES 

C 

50  00  10  K»l,9 
L-K-l 

C THE  FOLLOWING  VALUE  IS  8 FACTORIAL 

XF8»40320.0 

C 

C ....THE  FOLLOWING  CALCULATES  L FACTORIAL  

C 

XFL  • 1.0 
00  15  Ll-l.L 
Cl  • FL0ATIL1I 
15  XFL  * XFL*Cl 
C 

C ....THE  FOLLOWING  CALCULATES  (0-LI  FACTORIAL 

C 

L8-8-L 
XF8L  • 1.0 
00  20  L2-1.I.8 
C 2 • FL0ATIL2) 

20  XF8L  - XF8L4C2 
00  1 I • 1 , NX 

1 IF  IM.EQ.NPOINTIDI  GO  TO  2 
SIG»0.0 

WRITE  16,2101  M 
RETLRN 
C 

C SET  UP  INITIAL  PARAMETERS  

C 

2 JFIRST-REVU) 

JLAST  • KEV(l«ll-l 
C-l  ft-A I 22.0 
0-(B*AI/2.0 


SUM-0.0 

C 

C SUMMATION  OF  TERMS 

C 

UO  5 J-JFIKST , JLAST 

IF  IZ ( J ) .LF.O.O)  SUM  - SUM  ♦ WEIGHTfJI  ♦ <2.*FX1(D)  * 

1 (1.  -F  XI  (01  I*  (XF8/(XF0L*XFU)  * F X2  ( 0 ) **L  ♦ ( l . -FX2  ( 0 ))  • *(  0-1  I *0  > 

5  IF  (2UI.CT.0.0)  SUM  - SUM  ♦ WEIGHT(J)  ♦ ( ( 2 . * ( F X 1 < l < J I *C  *0  I • 

A (l.-FXK7(  J)*C*0))l  * (XF8/(XF0L»XFL ) ) ♦ ( FX2  ( Z ( J ) ♦OO  ) **L  • 

B (1.-FX2I Z( J)*C*D)  ) (8-L)  * (Z(J)*C«0))  ♦ ( 2 . *F  X 1 1 -Z  ( J ) *C  ♦O  ) 

C • (U-FXK-ZI JI*C*0)l  * I XF0/ (XF0L*XFL ))  • FX2 ( -! f J ) *C »D ) 

0 **L  * ( 1.-FX2C-Z ( J) ♦C+D) ) ( 8-L  l ♦ (-Z(JMOOII)) 


TWOPl-6. 283185 
SIG  - T WOP  1 • SUM 
10  WRITE  (6,202)  l, SIG 
GO  TO  70 


C 

C THIS  PORTION  PERFORMS  A GAUSSIAN  INTEGRATION  OF 

C Of  ANY  POLYNOMIAL  UP  TC  THE  FIFTH  OROER 

C 


60  00  6 I - 1 # NX 

6 IF  IM.EQ.NPOINTI I))  GO  TO  7 
SIG-0.0 

WRITE  (6,210)  M 

210  FORMAT  ( 1H0,10X, I2.2X.25HI S AN  INVALID  VALUE  FOR  M) 

RETURN 

C 

C SET  UP  INITIAL  PARAMETERS  ....... 

C 

7 JF I R S T-KE  Y ( I ) 

JLAST  - KEYII ♦11—1 
C ■ ( B-A ) /2 • 0 
0-(B>A)/2.0 

SUM  *0. 0 
C 

C SUMMATION  OF  TERMS 

C 

DO  8 J-JFIRST, JLAST 

IF  (Z(  JI.LE.O.O)  SUM  - SUM  ♦ WEIGMTCJ)  ♦ FXKO)  • 0 

0 IF  (Z(J).GT.O.O)  SUM  - SUM  ♦ WEIGMT(J)  * ((FXl(Z(J)*C  ♦ 01  • 

1 (Z(J)  * C ♦ Dll  ♦ ( F XI (~Z ( J ) • C ♦ D)  * (-Z(J)  ♦ C ♦ 0)1) 

SIG  • OSUM 
WRITE  (6,205)  SIG 
200  FORMAT  ( INI ,10X,I3A6) 

202  FORMAT  ( 1H0,10X,21HTHE  CROSS  SECTION  FOR , 2X 1 1 , 2X , 20HL-SHELL  VACANC 
lies  IS,2X,F10.6) 

205  FORMAT  ( IMO, 10X.29HTHE  SUM  OF  THE  INTEGRATION  IS.F10.6) 

70  RETURN 
END 


flBFTC  0B 

C THIS  SUBROUTINE  DESCRIBES  THE  K-SHELL  PROBABILITY  OF  IONIZATION 

FUNCTION  F X 1 (X) 
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COMPUN  Sl.S2.S3.SA.S5.S6 
COMPON  Tl,T2.T3,TA,T5,T6 

FX1"S6»X»»5  ♦ S5*X»*A  » SA»X**3  ♦ S3»X*»2  ♦ S2»X  ♦ SI 

RETURN 

END 


MBFTC  CC 

C THIS  SUBROUTINE  OESCRIBES  THE  l-SHELl  PROBABILITY  OF  IONIZATION 

FUNCTION  FX2  C X I 
COUPON  S1,S2.S3,S*,S5.SS 
COUPON  T1,T2,T3,T4,T5,T6 

FX2*T6»X»»5  ♦ T5»X»*A  ♦ T**X»»3  ♦ T3»X«*2  ♦ T2»X  ♦ Tl 

RETURN 

END 

THE  FIRST  LOCATION  NOT  USEO  BY  THIS  PROGRAM  IS  00660. 

THE  FIRST  LOCATION  NOT  USED  BY  THIS  PROGRAM  IS  01217. 

THE  FIRST  LOCATION  NOT  USED  BY  THIS  PROGRAM  IS  00113. 

THE  FIRST  LOCATION  NOT  USED  BY  THIS  PROGRAM  IS  00113. 

ICATA 


SC  A 

CALCULATION 

HE 

ON 

AL 

Vl/VL  - 3.0 

THE 

INPUT 

FUNCTION  FOR 

FXl  IS 

-0.  *X**5  ♦ 

-0.000641*X**4  ♦ 

0.C05047*X**3  ♦ 

-0. 

008975*X**2  ♦ 

- 

0.012925*X  ♦ 0.036866 

THE 

INPUT 

FUNCTION 

1 FOR 

FX2  IS 

-0.  *X**5  ♦ 

-0.  *X**4  ♦ 

0.000475*X**3  ♦ 

-o. 

003348* X**2  ♦ 

- 

0.000289*X  ♦ 0.029544 

THE 

INPUT 

values  FOR  A 

tB# 

ANO  M ARE » 0.200000 

4.000000  ANO 

15 

THE 

CROSS 

SECTION 

FOR 

0 

L-SHELL 

VACANCIES 

IS 

0.419104 

THE 

CROSS 

section 

FOR 

1 

L-SMELL 

VACANCIES 

IS 

0.072103 

THE 

CROSS 

SECTION 

FOR 

2 

l -SHELL 

VACANCIES 

IS 

0.006051 

THE 

CROSS 

SECT'ON 

FOR 

3 

L-SHELL 

VACANCIES 

IS 

0.000307 

the 

CROSS 

SECTION 

FOR 

4 

L-SHELL 

VACANCIES 

IS 

0.000010 

the 

CROSS 

SECTION 

FOR 

5 

L-SHELL 

VACANCIES 

IS 

0.000000 

The 

crc*„ 

SECTION 

FOR 

6 

L-SHELL 

VACANCIES 

IS 

0.000000 

the 

CROSS 

SECTION 

FOR 

7 

L-SHELL 

VACANCIES 

IS 

0.000000 

the 

CROSS 

SECTION 

FOR 

8 

L-SHELL 

VACANCIES 

IS 

c'.oooooo 

76 


APPENDIX  B 
SENSITIVITY 

The  measurements  of  the  A1  and  Si  Ka  x-ray  spectra  can  be  used  to 
determine  the  dependence  of  the  sensitivity  of  x-ray  detection  with  a 
crystal  spectrometer  on  the  atomic  number  and  energy  of  the  projectile. 
Sensitivity  is  nf  utmost  importance  in  analytical  applications  of  these 
measurement  procedures  and  in  measuring  the  Ka  x-ray  spectra  to  study 
the  satellite  dependence  on  such  factors  as  the  chemical  environment 
and  the  atomic  number  of  the  target  atom.  In  analytical  applications, 

Ka  x-ray  spectral  information  can  be  used  in  quantitative  analysis,  in 
the  determination  of  chemical  composition,  and  in  measuring  the 
thicknesses  of  surface  layers. ^ 

The  curves  in  Fig.  17  give  a relative  measure  of  the  dependence  of 
the  sensitivity  of  x-ray  detection  on  the  atomic  number  and  energy  of 
the  projectile  and  illustrate  the  importance  of  proper  selection  of 
projectiles  and  projectile  energies  for  maximum  sensitivity.  The  ratios 
of  the  total  A1  and  Si  Ka  x-ray  peak  areas  to  the  total  background  areas 
(beneath  the  x-ray  peaks)  in  spectra  obtained  with  He,  C,  0,  and  Ne  ions 
ranging  in  energy  from  1 to  10  MeV/amu  are  shown  in  Fig.  17.  The  curves 
for  aluminum  and  silicon  are  basically  similar.  At  the  lower  energies, 
all  the  curves  exhibit  peak  to  background  (P/B)  ratios  between  5 and  20 
which  rapidly  drop  as  the  projectile  energy  increases.  Between  2.5  and 
4.0  MeV/amu,  the  P/B  ratios  abruptly  change  to  lower  rates  of  decrease 
with  increasing  projectile  energy.  The  rapid  drop  in  the  P/B  ratios 


\ 


• 


FIG.  17.  The  dependence  of  the  peak  to  background  (P/B)  ratio  for  A.1 
and  Si  Ka  x-rays  on  projectile  energy  and  atomic  number. 


e^/amu ) 
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occurs  in  the  energy  region  which  closely  corresponds  to  the  center-of- 
mass  projectile  energies  necessary  to  overcome  the  nuclear  Coulomb 
barriers  in  aluminum  and  silicon  (shown  by  the  vertical  arrows  in  Fig. 
17).  This  suggests  that  the  majority  of  the  background  measured  at  the 
higher  projectile  energies  is  the  result  of  nuclear  reactions.  The 
error  bars  shown  on  the  helium  ion  curves  in  Fig.  17  illustrate  the 
amount  of  variation  of  the  P/B  ratios  observed  in  several  different 
runs.  These  variations  can  be  caused  by  factors  such  as  beam  spot  size 
and  soller  slit  and  crystal  alignment,  and  may  well  account  for  the 
inconsistencies  of  the  carbon  on  aluminum,  helium  on  aluminum,  and 
oxygen  on  silicon  monoxide  curves.  The  oxygen  ion  data  were  obtained 
from  previous  measurements  made  by  Watson  et 

From  Fig.  17,  it  appears  that  overall  the  sensitivity  is  best  with 
neon  ions,  however,  the  sensitivities  for  Si  Ka  x-rays  are  not 
significantly  different  for  any  of  the  projectiles  below  approximately 
2 MeV/amu.  In  all  cases,  the  P/B  ratios  reach  their  highest  values  at 
the  lowest  energies  available  in  this  work  (0.5  to  1.5  MeV/amu).  This 
is  also  the  energy  region  at  which  the  L-shell  ionization  reaches  a 
maximum  (see  Fig.  16). 
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